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ABSTRACT 

Aims. Disk-averaged infrared spectra of Neptune between 1.8 and 13 nm, obtained by the AKARI Infrared Camera 
(IRC) in May 2007, have been analysed to (a) determine the globally-averaged stratospheric temperature structure; (b) 
derive the abundances of stratospheric hydrocarbons; and (c) detect fluorescent emission from CO at 4.7 fim. 
Methods. Mid-infrared spectra (SGI and SG2 channels of AKARI/IRC), with spectral resolutions of 47 and 34 re- 
spectively, were modelled using a line-by-line radiative transfer code to determine the temperature structure between 
1-1000 /ibar and the abundances of CII4, CII3D and higher-order hydrocarbons. A full non-LTE radiative model was 
then used to determine the best fitting CO profile to reproduce the fiuorescent emission observed at 4.7 /im in the NG 
channel (with a spectral resolution of 135). 

Results. The globally-averaged stratospheric temperature structure is quasi-isothermal between 1-1000 fihar, which sug- 
gests little variation in global stratospheric conditions since studies by the Infrared Space Observatory a decade earlier. 
The derived CII4 mole fraction of (9.0 ± 3.0) x 10~* at 50 mbar, decreasing to (0.9 ± 0.3) x 10~* at 1 /ibar, is larger 
than that expected if the tropopause at 56 K acts as an efficient cold trap, but consistent with the hypothesis that 
CH4 leaking through the warm south polar tropopause (62-66 K) is globally redistributed by stratospheric motion. The 
ratio of D/H in CH4 of 3.0 ± 1.0 x 10~* supports the conclusion that Neptune is enriched in deuterium relative to the 
other giant planets. We determine a mole fraction of ethane of (8.5 ±2.1) x 10~^ at 0.3 mbar, consistent with previous 
studies, and a mole fraction of ethylene of 5.0^2 ^ 10~^ at 2.8 /ibar. An emission peak at 4.7 /im is interpreted as a 
fluorescent emission of CO, and requires a vertical distribution with both external and internal sources of CO. Finally, 
comparisons to previous L-band studies indicate significant variability of Neptune's flux densities in the 3.5-4.1 /tm 
range, related to changes in solar energy deposition. 

Key words. Planets and satellites: individual: Neptune — Infrared: solar system 
1. Introduction order hydrocarbons (Fig. [1]). Neptune's bulk C/H ratio 



is known to be enr iched over the solar abundance of 
Disk-integrated spectra of Neptune between 1.8 and 13 Grcvcssc ct al. (2007^ by 54±15 (from the measurements of 
microns were recorded on May 13, 2007, using both the Brines ct al. 1995;). This is consi s tent with formation mod - 
prism and the grism modes of the Infrared Camera (IRC, (gg^^ g_g iHersant et al. 2004 iQwen fc Encrenaz 200fil^ . 
l Onaka et al. 20071 ) on board ISAS/JAXA's AKARI m- ^hich predict that Neptune received a larger proportion 
frared astronomy satellite (|Murakami et al. 20071 ), which j^g^vy elements (relative to the Ha-He envelopes) com- 
launched on February 21, 2006 (UT). The broad p^^d to the other giant planets. AKARI spectroscopy 
wavelength-coverage provided by AKARI permits simul- ^oth CH4 and the higher-order hydrocarbons will 
taneous observation of near-IR reflectance of sunlight and ^^^^ ^o determine Neptune's stratospheric tempera- 
mid-IR thermal emission. These data are used to constram ^^^^ composition. The cold atmospheric tempera- 
models of Neptune's atmospheric structure, composition ^^^.^^ ^^^y^ an extremely low-power thermal emission spec- 
and aerosols in the troposphere and stratosphere. trum which requires high sensitivity to measure accurately. 
Interior models indicate that Neptune's atmospheric en- and AKARI results will be compared to previous deriva- 
velope of H2 and He represents only a thin shell com- tions of Neptune's thermal struc ture from the Voy ager 
pared to the radius of the planet (e.g. |Guillot 199^), but radio science investigation (RSS. iLindal et al. 1990fl . in- 



the collision-induced continuum of Ha-He is the dom- feared spectrometer (IRIS, IConrath et al. 19981) and the 

inant source of opacity in the mid-IR spectrum (6-13 Infrared Space Observatory riSO.lFeuchtgruber et al. 19991: 

/im). In addition to continuum emission, the 6-13 fim iFouchet et al. 2003HB"urgdorf et al. 2003D . 
spectrum exhibits substantial structure due to the chem- 
istry of stratospheric methane (CH4), CH3D and higher- 
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Previous studies have shown that the stratospheric 
abundance of CH4 is larg er than th e saturated abun- 
dance at the tropopause (B aines fc Han imel 1994), which 
should serve as an efficient 'cold trap' for CH4, restrict- 
ing it to the troposphere. Vigorous vertical advection of 
both gaseous CH4 and CH4 ice (which subsequently sub- 
limates) was invoked to explain the stratospheric abun- 
dance. Such rapid convection is consistent with observa- 
tions of Neptune's strong internal heat flux, t he disequi- 
librium of para- hydrogen (jConrath et al. 19981 ) . the high 
levels of tropospher ic CO determined from the sub-mm 
(jHesman et al. 20071 ) and with the visible observations of 
thick, high clouds casting shadows on deeper levels in 
Voyager 2 images. However, convective cloud activity is 
typically restr icted to Neptune's mid -latitudes and is time- 
variable (e.g. ISromovskv et al. 20011) . and it remains un- 
clear whether the strength of the upwelling is sufficient to 
maintain the stratospheric CH4, which is continually de- 
stroyed by p hotolysis. As an alt ernative, ground-based mid- 
IR imaging (lOrton et al. 200?! ) demonstrated the presence 
of a hot south polar region, warm enough to permit signifi- 
cant CH4 escape into the stratosphere without the need for 
strong convection. Ultimately, it is likely that both mech- 
anisms contribute, but AKARI derivations of the global 
CH4 abundance will be used to reassess the requirement 
for strong convection of methane-laden air into the strato- 
sphere. 

Once in the stratosphere, CH4 is photolysed to form 
the higher hydrocarbons. Acetylene (C2H2 at 13.7 fim) and 
ethane (C2H6 at 12.2 /im) were first de tected in ground- 
based spectroscopy (lOrton et al. 19871) and confirmed 
by Voyager/IRIS spectra ( Bezard et al. 199TI) . Ethylene 
(C2H4 at 10.5 n m) was detected in ISO spectra by 
ISchulz et al. (1993 ). and both methylacetylene (C3H4 at 
15.8 /zm) and diacetylene (C4II2 at 15.9 u.m) were recentl y 
detected in Spitzer/IRS spectra ()Meadows et al. 2008[ ). 
These hydrocarbons may diffuse downwards and freeze out 
as haze layers, and sunlight absorption on these hydrocar- 
bon hazes may be partially responsible for the warm tem- 
peratures of the lower stratosphere. 

The near-IR region of the AKARI spectrum (1.8-5.0 
/im) is dominated by the reflection of sunlight from cloud 
and haze layers, with broad absorptions due to atmospheric 
CH4. Clouds of CII4 are expected to condense at the 1-2 
bar level, overlying much deeper clouds of NH4SH and II2O, 
but only two of the uppermost layers have been observed: a 
thin cloud near 1.5 bar and an optically thick one with cloud 
tops near 3.8 bar (from ground-b ased H2 quadrupole lines 
and CH4 reflectari ce spectroscopv iBaines fc: Hammel 19941 
iBaines et al. 1995D . The near-IR spectrum will be used to 
determine the vertical distribution of CO from a fluorescent 
peak detected at 4.7 /im. 

The AKARI/IRC data reduction methods for both the 
near-IR and mid-IR spectra are described in Section [5] In 
Section[2]we use the mid-IR spectra to determine the global 
stratospheric temperatures, the abundances of CII4 and 
CH3D and estimates for the hydrocarbon mole fractions. 
The temperature structure is used in Section H] to model 
the fluorescent emission of CO. The implications of each of 
these new compositional measurements will be discussed in 
Section [S] 



2. Observations 

Fig. [T] shows the near and mid-IR spectra of Neptune on 
May 13, 2007, with many of the emission and absorption 
features described above. The emission features of CH4 and 
C2II6 are particularly prominent in the mid-IR spectra, 
with CII4 absorptions visible in the reflected sunlight com- 
ponent of the near-IR. In addition, we label peaks which 
may include some Hj^ emission (although these cannot be 
quantified properly in the low-resolution AKARI spectra) 
and a fluorescent emission due to CO at 4.7 /xm (see Section 

These disk-averaged spectra were obtained in four mod- 
ules (TablelH): SGI (5.5-8.3 nm) and SG2 (7.4-13 fim) chan- 
nels in the mid-IR, and prism (NP, 1.7-5.5 /zm) and grism 
modes (NG, 2.5-5.0 /im) in the near-IR. The spectral resolv- 
ing power was 135 in the interval 2.5 - 5 //m and between 20 
and 50 at other wavelengths. The two grisms on the mid-IR 
channel observe simultaneously, whereas either the NP or 
NG modes must be selected on the near-IR channel. The 
lower-resolution near-IR prism spectrum is not shown in 
Fig. [T] except for the 5.0-5.5 /im region, but the measured 
flux was found to be consistent between the NG and NP 
observations except between 4.1 and 4.7 /im, where the val- 
ues from the prism are systematically some 20% higher than 
those obtained with the grism. This may simply be due to 
an error in the flux calibration of the slitless spectrum or 
it may imply temporal or spatial variations in Neptune's 
atmosphere (Section [53]) • Long-term variations by a factor 
of two in the L band have been observed already in the past 
(|Feuchtgruber fc Encrenaz 20031) . 

Raw data were reduced using version 2.0 of the IRC 
spectral reduction pipeline, which includes masking soft- 
ware to avoid saturation in the mid-IR at bright hydro- 
carbon emission features. Each pointing consisted of eight 
exposures (nine for the grism in the near IR) , and any sub- 
frames with systematic noise patterns were removed from 
the analysis. As Triton was included in the field of view at 
a distance of only 6 arcsec in the cross-dispersion direction, 
we had to subtract its contribution to the near-infrared 
spectra. For the prism and grism modes, Triton's contri- 
bution was obtained from off-source spectra, extracted at 
the same distance from the moon as Neptune but on the 
other side, and then subtracted. The contribution of Triton 
to the prism mode was found to be much larger than to 
the grism mode, because in the slitless spectrum Triton 
stood "at longer wavelengths" relative to Neptune in the 
dispersion direction. As a consequence, the subtraction of 
Triton from the NG spectrum would have only served to in- 
crease the spectral noise, so was omitted in this reduction. 
NG shows good agree ment with ISO at 2.7 fim (based on 
lEncrenaz et al. 2000bD . 

The contribution due to Triton was deemed negligible 
in the mid-IR, so it has not been subtracted from the spec- 
tra. Spectra were extracted from a 3-pixel wide column for 
SGI and 5 pixel column for SG2. The short-wavelength end 
of SGI (shortward of 5.8 /im) demonstrated a rising flux 
which was deemed spurious (because the signal-to-noise ra- 
tio was rather low in this region), hence the 5.5-5.8 /im 
region is omitted from Fig. [T] As the nominal range of the 
NP spectrum ends at 5.2 /xm, there is a discontinuity in the 
spectrum between five and six microns. 
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2.1. Field of View 



3.1. Reference Neptune Atmosphere 



AKARI does not spatially-resolve Neptune's disc, which 
subtended 9.45 x 10~^^ steradians at the time of the obser- 
vations. With the exception of the NG grisni mode (where 
we used a sUt 1 arcmin wide) , all spectra were obtained in a 
slitless mode (lO'xlO') so that Neptune was fully included 
in the FOV. The fact that Neptune is an extended source 
(2.3 arcsec compared to a 1.5-arcsec pixel size for the NIR 
camera) was taken into account during the reduction. For 
the near-IR spectra of Neptune we used 3 pixels along the 
cross-dispersion axis and checked that there was negligi- 
ble signal outside the extraction rectangle. 3 pixels are the 
default (this is the FWHM of the inflight point spread func- 
tion in the near-IR) , and using a wider column would have 
further contaminated the spectrum with flux from Triton. 
In the mid-IR, Neptune's diameter is even more negligible 
compared with the instrumental PSF than in the near-IR. 



2.2. Wavelengtli Calibration 

The pipeline determines the location of various sources and 
uses a linear relation between pixel location and wave- 
length. The slitless nature of the AKARI/IRC, in addi- 
tion to observing an extended source and the problem of 
saturation in longer exposures which can mask the precise 
location of a source, can lead to difficulties when determin- 
ing the wavelength calibration. Furthermore, as spectra can 
only be extracted from the array on an integer-pixel basis, 
this can introduce systematic pixel errors as large as ±0.03 
fim. The positional uncertainty for wavelength calibration 
in the NP spectrum is also of the order of a single pixel, 
leading to a wavelength calibration offset of approximately 
0.05 /im. 



2.3. Instrument Function 

Information on the spectral lineshape was not available at 
the time of writing, so we use the assumption of a Gaussian 
with the required spectral resolution, R — A/AA, where 
AA is the smallest resolving element, which should be a 
constant in each module to a first approximation. A com- 
parison between spectra computed with a constant R and 
a constant AA indicated that this approximation was valid. 
The specified resolutions are given in Table [T] 



3. Analysis of Mid-IR Thermal Emission: 5.8-13 fim 

Spectra from the AKARI/IRC SGI and SG2 mid- 
IR channels were analysed using an updated ver- 
sion of the line-by-line radi ative transfer algorithm, 
RADTRANS (jlrwin et al. 19971 ). The disk-averaged trans- 
mission was calcu lated via an exponential integral method 
(jGoodv fc Yung 19891. Because the longest wavelength ac- 
cessible to AKARI is 13 /im, we are unable to use the 
smooth collision-induced II2 and He continuum longwards 
of 15 fj,m to perform a full derivation of the tropospheric 
temperature structure beneath 0.1 bar. However, AKARI 
spectra are sensitive to the stratospheric temperatures in 
the 1-1000 /ibar region via the emissions of the CH4 and 
the higher-order hydrocarbons. The following sections de- 
scribe the determination of stratospheric temperature and 
composition from the disk-averaged spectra. 



As AKARI spectra offer no constraint on the at- 
mosphe ric temperatures beneath the tropopause, we 
follow iMeadows et al. [^OS) and use the Neptune 



T{p) structure from iMoses et al. f20"05l) as a start- 
ing point for fitting temperatures above the 1-mbar 
level. The initial profile was formed from a combina- 
tion of results from Voyager radio science inv estiga- 
tion in the troposphere (RSS, iLindal et al. 19901 ). with 
updates in the stratosphere from the Infrared Space 
Observatory (iBez ard et a l. 19981 iFe uchtgr uber et al. 19991: 
iFouchet et al. 200 3: Burg dorf et al. 2 003) an d from 17-24 
^m ground-based spectra ( Orton et al. 198711 19901 ). The 
resulting stratospheric profile became isothermal at 168 K 
above the .1-mbar level, but was forced to 250 ± 50 K 
at 0.1 /zbar (jBezard et al. 19981 and references therein) for 
consistency with stellar occultation data. The profile was 
regridded onto a pressure scale between 10 bar and 0.1 /ibar 
in 100 intervals equally spaced in log(p). The gravitational 
acceleration at 1 bar is 11.2 ms^^ when we take a global 
average (from equator to pole, where g is strongest), and 
the forward model accounts for the variation of Neptune's 
gravity with latitude. The H2 mole f raction was set to 0.8 3 
with a helium mole fraction of 0.15 (jConrath et al. 19931 ). 

The reference temperature profile used in this anal- 
ysis had been used as input for a one-dimensional 
model o f Neptune's strato spheric photochemistry (Model 
A from iMoses et al. 20051) which provided the hydrocar- 
bon abundance profiles for scaling analysis in this data. 
The deep CH4 abunda nce was set to 2.2% following 
iBaines fc Hammel fl994l) who determined this value using 
a combination of the Voyager/RSS temperature profile and 
ground-based measurements of H2 quadrupole lines. The 
vertical CH4 reference profile follows a saturated vapour 
pressure profile with 100% humidity in the tropopause re- 
gion, but then increased back to a mole fraction of 2.2% 
when the stratosphere became warm enough. Above this 
level, the CH4 abunda nce decreased accor ding to the pho- 
tochemical modeling of IMoses et al. (2005[ ). Vertical hydro- 
carbon distribution s were included following the analysis of 
IMoses et al. ('20051 ) . and scaled to reproduce the observed 
spectra. We hold ^^C/^'^C constant at the terrestrial value 
of 89. Because AKARI spectra are insensitive to the S(0) 
and S(l) lines in the far-IR, the ratio of ortho-H2 to para- 
H2 was set to thermochemical equilibrium at every level 
of the atmosphere. We used this reference atmosphere to 
calculate the contribution function (the product of the rate 
of change of the transmission with altitude with the black 
body function), which shows the level of emission for each 
of the spectral regions covered by the AKARI spectrum in 
the 0.3-3000 /ibar range (Fig. The sources of spectro- 
scopic line data are the same as those used in Jupiter and 
Saturn analyses ([Fletcher et al. 200"9allbl ). 



3.2. Methane Emission at 7. 7 ^m 

The spectral resolution of the AKARI mid-IR data does 
not permit characterisation of the stratospheric tempera- 
ture independently of composition. Instead, we must rely 
on a trade-off analysis between two strongly anti-correlated 
variables - the stratospheric temperature and the methane 
mole fraction, which both determine the strength of the 
emission of the 1^4 CH4 band at 1304 cm~^. Such a correla- 
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tion leads to large error estimates on the derived quan- 
tities. Furthermore, the vertical variation of CH4 and 
the other hydroca rbons must be fixed to the models of 
iMoses et al. f2005l) . because we have insufficient spectral 
resolution to determine the vertical variability of these 
species. However, as the following analysis shows, the in- 
dependent measure ments of AKARI co nfirm the validity of 
the profiles used bv lMoses et al. C2005I ). 

Our approach to the derivation of T{p) was a fam- 
ily of single- variable parameterisations, each one building 
upon the last (Fig. [3]) in a bottom-up direction (i.e. start- 
ing with the 1-mbar level where the temperature is best 
known from previous determinations, and working our way 
to 1 ^hai and above). This approach was guided by the 
multi-lobed nature of the contribution function in Fig. [5J 
which shows sensitivity throughout the 1-1000 fihav range. 
We attempt to pre serve the smooth nature of T(p) inferred 
from ISO spectra (iBezard et al. 1991 iFouchet et al. 2003t 
iBurgdorf et al. 2003[) . but we caution that the presence of 
real oscillations in the temperature profile (due to strato- 
spheric wave activity, for example) would not be detected 
via this simple parameterisation. At the same time, we 
scaled the CH4 abundance above the tropopause and min- 
imised the fit to the data, x^j over the 7.1-8.3 /im (1200- 
1400 cm~^) spectral range. The four temperature profile 
parameterisations were as follows; 

1. Isotherm: Starting from a smooth extrapolation be- 
tween the 1 mbar level (approxi mately 136 K) and 0.1 
/xbar (approximately 250 K from lYelle et al. 1993[ ). we 
generate a fully isothermal stratospheric temperature 
at a range of test values between 130 and 250 K (Fig. 
[31[a)). The best-fitting solution required 160-165 K and 
a CH4 mole fraction of 1.0 — 3.5 x 10~^ at 50 mbar. 

2. Lapse Rate: Starting from the best- fitting isothermal 
structure, we add small positive and negative dT /dp 
gradients as shown in Fig. [31[b) , and smoothing the ver- 
tical structure so that there is some small perturbation 
at pressures greater than 1 mbar. The best-fit solution 
required a warmer temperature profile increasing from 
166 K at 100 /xbar to 170 K at 1 /ibar, and a smaller 
CH4 abundance of (1.0 ± 0.5) x lO'^ at 50 mbar. 

3. Smoothing the Profile: To further test the validity of a 
quasi-isothermal temperature structure, we use the fam- 
ily of T{p) profiles in Fig. ^c) to interpolate smoothly 
over the 1-1000 /ibar region. However, we found that 
the spectral fits were considerably worsened by the high 
temperatures at the //bar level, so these models were 
omitted from the subsequent analysis. 

4. Extrapolation to 0.1 /xbar: The model from step two 
reproduce the spectrum everywhere except for the cen- 
tre of the CH4 Q-branch, which probes the highest al- 
titudes (Fig. [2])- Preserving the quasi- isotherm in the 
10-1000 /ibar region, the profiles in Fig. |3Kd) were used 
to test our sensitivity to the high 0.1 /xbar temperatures 
(approximately 250 K) determined from stellar occupa- 
tions. Satisfactory solutions were found for a transition 
from the quasi-isotherm to the warm temperatures at 
altitudes between 0.1-5.0 /ibar, with a CH4 mole frac- 
tion of (6.0 ± 3.0) X 10-"* at 50 mbar. 

An example of the analysis is shown in Fig. 01 which 
shows a two dimensional x^ surface in the centre, the spec- 
tral effect of varying the CH4 mole fraction on the left and 
the temperature profile on the right. The two correlated 



parameters can be marginally separated because of the dif- 
ferent morphologies of their effects on the spectrum. The 
best-fitting temperature structure and the error range (tak- 
ing the anticorrelation with CH4 into account) are shown 
in Fig. [5l compared to T{p) determinations from previous 
studies. The warm and cool extrema of the error range were 
used in the error analysis for composition in the following 
sections. The best-fitting temperature profile provides an 
optimum CH4 mole fraction of (9.0±3.0) x 10"'' at 50 mbar, 
decreasing to (0.9 ± 0.3) x lO""* at 1 /xbar. However, using 
the cool and warm extrema of the temperature profile, the 
CH4 mole fraction could vary between 6.0 — 18.0 x lO"''. 
The temperature and CH4 determinations will be discussed 
in Section [S] 

3.3. CH3D at 8.6 fim 

Fig. ini shows the effect on the variation of the D/H ratio in 
CH4 on the AKARI spectrum. The vq vibrational band of 
CH3D at 1161 cm~' is seen entirely in emission, afi'ecting 
the shape of the spectrum on the short-wavelength side of 
the SG2 spectrum. Unfortunately, the absolute calibration 
(both in terms of flux, and the wavelength calibration) has 
the largest degree of uncertainty at the end of the channel: 
spectral reductions optimized for the two prominent mid-IR 
features, C2H6 and CH4, provided differing estimates of the 
flux in this region, so the optimization for CH4 had to be 
used. Furthermore, as no spectral lines are resolved, this can 
provide no more than an order of magnitude estimate of the 
CH3D abundance. Nevertheless, we scaled the CH4 profile 
by a CH3D/CH4 ratio between 2.0 x lO"'' and 6.0 x lO"'' 
in steps of 0.1 x lO"**, using the cool, warm and best-fitting 
T{p) structures from the previous section. 

Minimising x^ over the 1100-1200 cm~^ spectral range, 
we find that D/H ratios in CH4 in the range (3.0 ± 1.0) x 
lO"** can flt the spectrum, consistent with the value of 
3.25 x 10"'' that we would ex pect from the ISO-derive d 
D/H ratio in H2 of 6.5 x 10"^ (|Feuchtgruber et al. 1999t ). 
using a fractiona tion factor of / = 1.25 between deuterium 
in H2 and CH4 (jFedev fc Prinn 1988D . This D/H ratio in 
CH4 is pertinent to the 0.2 mbar level, the peak of the 
CH3D contribution function in Fig. [2] 

3.4. Hydrocarbon Emission Features 

The second most prominent feature of the mid-IR AKARI 
spectrum is the ethane peak at 12.2 fim, which is fltted in 
Fig. [71 Using the three temperature proflles, we scale the 
C2H6 profile in steps of 0.1 between 0.0 and 2.0, to derive 
a best fit of 0.8 ± 0.2 times the a priori profile. The error 
range was calculated using the three models representing 
the centre and the two extrema of the temperature pro- 
flles, yielding 0.6 for the warmest profile, 1.0 for the coldest 
profile, and 0.8 as the optimum profile. This yields a max- 
imum abundance of (2.5 ± 0.6) x 10~^ at 3.4 /ibar, and 
(8.5 ± 2.1) X lO"'^ at the peak of the ethane contribution 
function at 0.3 mbar in Fig. [21 

The long-wavelength end of the SG2 spectrum is sen- 
sitive to the abundance of acetylene (C2H2), where we see 
some discrepancy between the model and the data in Fig. [71 
Moses ct _al._i2QQ3) noted the difficulties associated with si- 
multaneously modelling the abundances of C2H6 and C2H2, 
suggesting that the model overpredicted the abundance of 
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acetylene. Indeed, a comparison of the model with data in 
the 760-780 cm~^ range implies the need for a scaling factor 
of 0.5-0.6 to reproduce the observations, suggesting a C2H2 
mole fraction of approximately 1.5 x 10^^ at 10 /ibar, or a 
peak abundance of 2.2 x 10"'' at 3.4 /xbar, though we stress 
that our sensitivity to C2H2 is limited to the wings of the 
emission band at 760-790 cm~^. 

Finally, the AKARI spectrum between 900-1100 cm^^ 
has a low signal-to- noise ratio (Fig.[T]), because the radiance 
is very weak at these wavenumbers, but shows a modulation 
at 952 cm~^ (10.5 fim, the 1^7 band) which could be caused 
by ethylene (C2H4). Fig.[5]indicates that the forward model 
allows us to partially reproduce the morphology of this re- 
gion of the spectrum. Using the thre e temperature pro - 
files, we determine a scale factor for the lMoses et al. I'2005D 
model profile of 1.4lQ g, equivalent to a peak abundance of 

l.Stli X lO"'^ at 1.3 /ibar, or 5.0t^ j x lO"'^ at 2.8 /xbar, 
the peak of the contribution function in Fig. [51 Though we 
must be cautious, because the spectral morphology is not 
precisely reproduced by the fo rward models. Furth ermore, 
the ISO spectrum presented bv lSchulz et al. Q999[ ) showed 
the C2H4 peak as a narrow single-pixel feature, whereas 
our modelled ethylene emission for AKARI is rather broad. 
Finally, the region of the spectrum between 980-1100 cm^-'^ 
shows excess emission which is not explained by our model, 
and it may be possible that higher-order hydrocarbons, 
whose spectral features are unresolved, are causing addi- 
tional emission in this range. 



4. Analysis of Near-IR Reflectance Spectra: 1.7-5.5 

Neptune's near-infrared spectrum in Fig.[T]is dominated by 
reflected sunlight from Neptune's cloud layers, in addition 
to absorption due to CH4. The peak at 4.7 fim is of particu- 
lar significance, as discussed in Section HTTJ At shorter wave- 
lengths, the AKARI near-IR spectrum shows considerable 
deviations from VLT/ISAAC spectra between 3.5-4.1 /zm 
obtained in August 2002 ( Feuchtgrubcr & Encr enaz 200311 . 
suggestive of meteorological variability discussed in Section 

4.1. CO Emission 

The emission peak at 4.7 fim cannot be reproduced by 
thermal or solar reflected LTE (local thermodynamic equi- 
librium) radiative transfer. Just as on Uranus, where 
a similar emission fe ature has been detected previously 
(|Encrenaz et al. 20041) . the feature can be interpreted as 
fluorescent (non-LTE) emission of CO. The population in 
the upper levels of the observed transition is by solar ab- 
sorption, mainly from the CO 1-0 and 2-0 bands. However, 
the contribution from the (1-0) band (resonant fluores- 
cence) is strongly self-absorbed, so the (2-1) band there- 
fore dominates the spectrum. This explains the wavelength 
shift in the observed spectrum compared to the (1-0) band 
centre. A full non-LTE radiative model has been adapted 
to Neptune with the following characteristics. The atmo- 
spheric structure was taken from the mid-IR determina- 
tions in Section [3l with a variable CO abundance vertical 
profile. Two CO vertical profiles were modeled to produce 
the synthetic spectra in Fig. [SI 



— Profile 1: CO limited to the stratosphere alone with a 
variable mixing ratio (constant with altitude), mixing 
ratio of CO/Il2=2.5 x 10^^, decreasing to zero below 
10 mbar. 

— Profile 2: Similar to the CO profile retriev ed from 
the sub-millimeter range (jHesman et al. 20071) ■ CO IS 
present in both the stratosphere and troposphere. 

The upper atmospheric abundance was mo deled using a n 
eddy diffusion coefficient of 2 x 10'' cm^s~^ (lAtreva 1992[ ). 
The synthetic spectra in Fig. [9] were calculated from a non- 
LTE radiative model including solar radiation absorption, 
self-absorption in the resonant fluorescent (1-0) band, and 
frequency redistribution from vibrational CO bands. For 
both profiles, we scaled the abundance of CO to obtain the 
best fit to the AKARI observations. 

The weighting function of the 2-1 fiuorescent emission 
is found to peak around 1 bar, with much smaller contribu- 
tions from the 1-0 band at higher altitudes (0.01-0.1 mbar), 
therefore giving sensitivity to the CO distribution in both 
the stratospheric and tropospheric parts of the atmosphere. 
A similar depth of emission was found for fluorescenc e 
of CO on Uranus, as discussed in lEncrenaz et al. f2004l ). 
Despite the fact that non-LTE effects are usually dominant 
only in the low-frequency coUisional regime of the upper 
atmosphere, the decrease in the non-LTE contribution is 
counter-balanced by the increase in the optical depth of 
CO, and the non-LTE emission in the 2-1 band is not af- 
fected by self-absorption like the 1-0 band. 

The results of the non-LTE radiative transfer model are 
shown in Fig. [HI A continuum has been added to the syn- 
thetic spectra to fit the AKARI spectrum outside the CO 
band. The dashed line shows the effects of restricting CO to 
the stratosphere - its absence from the troposphere means 
that the fiuorescent emission is insufficient to reproduce the 
peak in the AKARI spectrum. The best-fitting model had 
CO/Il2=2.5 X 10~^ in the stratosphere with a decrease in 
the abundance by a factor of four below 10 mbar. Only by 
including tropospheric CO do we begin to reproduce the 
emission feature, which would likely be improved via the 
addition of higher vibrations CO bands. However, given the 
spectral resolution of the observed spectra (and the uncer- 
tainties in the input atmospheric profiles discussed earlier), 
this was not deemed necessary at this stage. The conse- 
quences of tropospheric CO are briefiy discussed in Section 

m 

5. Discussion 

5.1. Globally- Averaged Temperatures 

Neptune's stratospheric temperature structure and com- 
position have been previously investigated by a number 
of authors, and in this section we show that the AKARI 
SGI and SG2 modules provide an independent verifica- 
tion of many of th ese results, particul arly of the photo- 
chemical models of iMoses et al. f 20051 ). Fig. [5] compares 
the derived stratospheric temperature to previous stud- 
ies, and demonstrates convergence of multiple models in 
the 0.1-100 /xbar range. The resulting temperatures are 
at the lower end of the 170-190 K temperature range de- 
rived from the Voya ger Ultrav iolet Spectrometer between 
0.3-50 /ibar fUVS. .Yelle et al. 1 9931. For pressures greater 
than 1 mbar the Voyager/RSS results stand out as be- 
ing cooler (minimum of 49 K, iLindal et al. 1990ll than 
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the others, where as subsequent analyses of remote sens- 
ms dat a from ISO ('Bczar d et al. 19 98VFouchct et al. 2003t 
[Burgd orf et al. 2 003) and the Voyager In frared Radiometer 
and Spectrometer (IRIS, IConrath et al. 19981 ) have sug- 
gested warmer tropopause temperatures of approximately 
56 K, dropping to 50-52 K at mid-latitudes. 

At pressures less than 1 mbar the AKARI quasi- 
isothermal temperature profile is remarkably consis- 
tent with that used in the photochemical analy- 
sis of Im OSes et al. ^20051 ). and warmer in the 0.1-1.0 
mbar region than the profile of iBezard et al. fl9"98l) . 
However, the AKARI profile is fully consistent with the 
ISO/SWS-derived temperature of 170 ± 3 K at 3 ^bar 
(jBezard et al. 19981 ). The SGI spectrum also requires that 
the transition from the quasi-isotherm to the warm tem- 
peratures at 0.1 fihai must occur at pressures lower than 
5.0 /ibar so that they do not adversely affect the quality of 
the spectral fits. The lack of sensitivity at these high alti- 
tudes is reflected in the large error bars above 1 /ibar. The 
close correspondence of all of the stratospheric temperature 
models suggests little variability in Neptune's global tem- 
perature structure in the decade between the ISO (1997) 
and AKARI (2007) observations. But disk-averaged spec- 
troscopy is insensitive to spatial variations in stratospheric 
temperatures across the planet (e.g. the localised south po- 
lar stratospheric emission discussed by Orton et al. 2007), 
and it is likely that the compositional and thermal results 
derived here are weighted towards the high-temperature 
south polar region. 

5.2. Atmospheric Composition 

Given the large anticorrelation between stratospheric tem- 
peratures and the CII4 abundance, it is perhaps surpris- 
ing that the temperatures match so closely when the 
methane mole fraction was varied as a free parameter. 
With saturation at the coolest tropopause temperatures 
of 56 K, we would expect a stratospheric mole fraction of 
2.2 X 10--*. The derived CH4 abundance of (9.0 ± 3.0) x 
10~^ at 50 mbar is four times as large as that expected 
from the cold-trap temperature. This result is consistent 
with the modelling of v isible CH4 absorption spectra by 
(|Baines fc Hammel 19941 ) . who found stratospheric abun- 
dance in the range 0.2 — 17.0 x 10^'*, with an optimum 
value of 3.5 x lO"** . Furthermore, the ranges of CH4 mole 
frac tions quoted bvlOrton et al. (19921 ) (1.9 - 26.1 x 10"^) 
and lOrton et al. f2007t) (7.5 - 15.0 x 10'"^) encompass the 
AKARI results at mbar pressures. Most importantly, the 
optimum value for CII4 is consistent within the errors with 
the stratospheric abu ndance of 7.5 x 10"*^ used in the photo- 
chemical modeling of iMoses et al. (2005|) . discussed below. 

Fig. (TU] compares the best-fitting abundances of CII4, 
C2H6 and C2II4 to the predict ed vertical profi les of the 
ID photochemical model of fjMoses et al. 20051 Fig. 32). 
The retrieved values have been placed at the peaks of 
the contribution functions. If we extrapolate the AKARI- 
derived CH4 mole fra ction to lower pres sures using the 
altitude-dependence of iMoses et al. f 20051 ). we find values 
of (0.9 ± 0.3) x 10""^ at 1 /xbar. This is at the upper limit 
of the 0.25 — 0.90 x 10""^ range of m ole fractions at 0.4- 
1.5 /ibar derived bv lYelle et al . (1993*) from Vo yager/UVS 
data (values extracted from|Mpses et al. 2005L and shown 
as a dotted box i n Fig. [TU]). Using a similar dataset, 
iBishoD et al. (19^^ found a 0.2 — 1.5 x 10~^ range between 



0.06-0.25 //bar, which is larger than the estimates from 
the AKARI-dcrived values with the model vertical profile 
(dashed box in Fig.fTUl). AK ARI results are more consistent 
with the lYelle et al. fl993|) stu dv. which improved on the 
analysis of iBishop et al. (19921 ) by using the full 126-166 
nm range of the Voyager/UVS data. So despite the non- 
uniqueness of the solution and the anti-correlation between 
the temperature and CII4 abundance, the AKARI thermal- 
IR results appear to be consistent with the UV-derived re- 
sults at higher altitudes, confirming the stratospheric super- 
saturation of CII4. This suggests that some form of vertical 
transport of CIl4-laden tropospheric air is required to sup- 
ply methane into the l ower stratosphere, as had been previ- 
ously postulated (e.g.. iLunine fc Hunten 19891) . This could 
be the result of vigorous convective activity on a global 
scale, but it is interesting to note Nept une's south pole is 
warm enougn (62-66 K at 100 mbar from lOrton et al. 20071 
compared with our nominal value of 56 K) to permit CH4 
mole fractions of 80 — 100 x 10^"* to be present in the 
south polar region, almost an order of magnitude larger 
than the global mean derived here. The warm south pole, 
along with elevated stratospheric CH4 , could dominate the 
disc-averaged spectrum. CH4 could be subsequently redis- 
tributed to other latitudes by stratospheric transport, but 
spatially-resolved spectroscopy of CII4 is required to distin- 
guish between the effects of localized convection versus po- 
lar enhancements. Disc-averaged AKARI observations are 
insufficient to distinguish between these two possible mech- 
anisms for the stratospheric enrichment in methane. 

We derive an estimate of the D/H ratio in CH4 of 
(3.0 ± 1.0) x 10~^ which is consistent with the estimate 
of (3.6 ± 0.5) X 10~"^ from g round-based spectroscopic 
observations of lOrton et al. (1 992)). The lower accuracy 
in our value reflects the uncertainties in the tempera- 
ture structure and CH4 abundance, along with the abso- 
lute flux uncertainties at the short-wavelength end of the 
AKARI SG2 channel. Nevertheless, this result is consis- 
tent with the value of 3.25^0 75 x 10"'' that we would 
expect from the ISO-derived D/H enrichment in H2 of 
X 10"^ (|Feuchtgruber et al. 1999l) . using a fraction- 
ation factor of / = 1.2 5 between deuterium in H2 and CH4 
(|Feglev k Prinn 19881 ). The AKARI analysis confirms the 
enrichment in Neptune's deuterium abundance relative to 
the other giant planets. 

The ethane mole fraction of (8.5 ±2.1) x 10" at the 
peak of the ethane contribution function at 0.3 mbar is 
smaller than the model predictions by approximately 20% 
(Fig. [To)) . Nevertheless, it is consistent with other infrared 
spectroscopic measurements of the C2H6 emission from het- 
erodyne techniques (1.9 x 10"^ with a factor of four uncer- 



tainty, lKostiuF|er^LT992), echelle spectroscopy (1.0 



+0.2 
-0.8 



10-6.lOrton et al. 1992[) and Voyager/IRIS {l.btH x 10" 



iBezard et al. 19911 ). Furthermore, the close correspondence 
between the photochemical model and the AKARI result in 
Fig. [TU] indi cates the consiste ncy between thermal-IR and 
UV results, l^e et al (1993!) used Voyager/UVS to show 
that the ethane abundance increased from approximately 
1 x 10"^ at 1 mba r to 3 x 10~^ at 1 / /bar, which was sub- 
sequently used by iMoses et al. (20051) as a constraint on 
the model curves for Fig. [10] Finally, this vertical profile 
was scaled in the AKARI analysis to obtain the C2H6 mole 
fraction presented here. 
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Because of the low spectral resolution of the AKARI 
spectra in the 10-/ini region, the abundance of ethylene 
at 2.8 /ibar (Fig. [TU)) should be considered with cau- 
tion. Our ethylene estimate of S.Olj'^ ^ 10~^ at 2.8 /ibar 
is approximately 40% larger than the model prediction. 
Extrapolatii ig to deeper press ures using the vertical de- 
pendence of lMoses et al. f2005[ ). we estimate the 200-/ibar 
mole fraction to be 1.4^Q'g x 10~^. This value is at the 
upper limit of the range of the C2H4 models tested by 
ISchulz et al. fl999 l) in their fitting of ISO photometry data 
for the same region. Indeed, Moses et al. (2005) 

point out that their ph otochemical mode l is more consistent 
with the UV results of lYelle et al. flQQSt ). which showed an 
increase from « 1.0 x 10~^ to ~ 5.0 x 10~^ between 100 
/ibar and 10 /ibar. Given that this vertical profile was scaled 
by 1.4^Q g to fit the AKARI spectra, we again find good 
agreement with the photochemical model and UV-derived 
results. 

An important test of photochemical models is their abil- 
ity to reproduce the observed ratios of the hydrocarbon 
species. The ethane/acetylene, ethane/ethylene and acety- 
lene/ethylene ratios should increase with distance from the 
Sun because of changing photolysis rates, but a combi- 
nation of th e AKARI results suggest that Model A of 
iMoses et al. (2005) overestimates the ratios by a factor of 2- 
3 at 0.2 mbar: ethane/acetylene is « 20 (model predicts 60), 
ethane/ethylene is « 680 (model predicts 1180) and acety- 
lene/ethylene is « 30 (model predicts 70). These difficulties 
encountered in reproducing thes e ratios in the photo chemi- 
cal models has been discussed bv lMoses et al. (20051 ). and it 
should be borne in mind that the estimates of hydrocarbon 
abundances depend strongly on the temperature structure, 
which is itself non-unique. High spectral resolution obser- 
vations of the hydrocarbon emission bands may help to 
resolve these open questions. 

At shorter wavelengths in the near-IR, the non-LTE 
modelling of the emission peak at 4.7 /im (Fig. [5]) demon- 
strated the need for the presence of CO in Neptune's tropo- 
sphere, independently con firming the interpreta. tion of sub- 
millimeter observations bv iHesman et al. f2007l ). Although 
the existence of CO in Neptune's stratosphere has been 
kn own since the first millimeter wavelength observations 
of iMarten et al. ('1993D . AKARI is the first instrument 
to observe the fluorescent emission of CO at 4.7 fj,m. 
Further more, the non-LTE modeling confirms the find- 
ings of E esman et al. (^2007t ) that an additional reservoir 
of CO is needed in the upper troposphere to fit the spec- 
tra. Thus Neptune's CO abundance is a combination of 
CO from the thermochemistry of carbon and oxygen in 
the deep atmosphere, advected to altitudes accessible to 
remote sensing, in addition to stratospheric CO from ex- 
ternal sources. Our CO distribution (CO/H2=2.5 x 10~^ 
in the stratosphere, decreasing by a factor of four below 
10 mbar) is entirely consistent with the best-fitting result 
of iHesman et al. (2007). The fiuorescent emission spectra 
are inconsist ent with the constan t vertical CO distribution 
proposed by iMarten et al. fl993h to fit the first measure- 
ments of C O emission in the m illimeter range. However, 
as noted bv IHesman et al. (20071 ). these early results were 
sensitive only to the emission cores, and not to tropospheric 
absorption further from the line centre, so that the constant 
vertical profile was a perfectly adequate assumption. The 



AKARI modeling confirms that such an assumption is no 
longer valid. 

CO is a disequilibrium species in the upper troposphere 
(it is not thermochemically stable at the temperatures and 
pressures observed), so must be supplied by convective 
mixing from the deeper atmosphere. On the other hand, 
stratospheric CO is probably the result of infalling exo- 
genic material, eithe r supplying CO dire ctly (a cometary 
impact hypothesis, iLellouch et al. 20051 ) or by supplying 
H2O to react with CH4 photochemical pr oducts (mete- 
oriti c bombardment. iRosenqvist et al. 1992D . As suggested 
by Hcs man et al. ('2007^. the determination of the CO pro- 
file alone is insufficient for discrimination between the dif- 
ferent external sources. Spatially resolved spectroscopy of 
CO and other 0-bearing species would provide invaluable 
clues to understand the stratospheric CO abundance. 

5.3. Meteorological Variations 



iFeuchtgruber fc E ncrenaz f2003l ) measured the spectrum 
of Neptune in August 2002 in the L-band (3.5-4.1 /im) with 
the ISAAC imaging spectrometer at the ANTU of the Very 
Large Telescope at the European Southern Observatory in 
Cerro Paranal. They found the observed flux to be stronger 
by a factor of three compared to measurements at the same 
wavelength recorded in 1997, five years earlier, and at- 
tributed the excess of infrared flux to an increased produc- 
tion of photochemical aerosols. This result suggested that 
Neptune's meteorology was actually linked to the variation 
of the solar incident flux. A similar conclusion was drawn 
from visible photometry of Neptune, which is correlated 
with the sol ar activity cycle via solar UV as a causative 
mechanism ([Lockwood fc Thompson 2002f) . In addition to 
these long-term trends, however, ther e were also iso lated in- 
frared outbursts, for example in 1976 (jPilcher 19771) . In con- 
trast, the AKARI measurements in the same spectral range 
produces flux densities 4 times smaller than those me a- 
sured 5 years earlier bv IFeuchtgruber fc Encrenaz f 20031 ). 

In Table [5] we compare the average flux density of 
Neptune in the L-band as measured by AKARI with the 
previous findings (estimated from spectra in the relevant 
publications), and with the number of sunspots (represen- 
tative of the solar variability for the incident flux). The 
correlation coefficient between fiux density and sunspot 
number for the latest three dates is 0.99. This is surpris- 
ingly high given the considerable uncertainties in both vari- 
ables, but supports the connection between L-band flux 
density (related to the generation of photochemical aerosols 
in the stratosphere) and solar i ncident flux discussed by 
IFeuchtgruber fc EncrenaT l'2003l) . 

6. Conclusions 

Infrared spectra from the AKARI Infrared Camera be- 
tween 1.8 and 13 /im have been analysed to determine 
the stratospheric temperature structure and composition. 
The 7.7-/im CH4 emission band was used to simultane- 
ously derive a quasi-isothermal temperature structure (con- 
sistent with previous studies in the 1-1000 /ibar range) and 
a CH4 stratospheric mole fraction of (9.0 ± 3.0) x 10^^ 
at 50 mbar, decreasing to (0.9 ± 0.3) x 10"'' at 1 /ibar. 
The close agreement in the globally averaged tempera- 
ture structure derived from ISO-observations in 1997 and 
AKARI-observations in 2007 suggest negligible changes to 
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Neptune's global stratospheric temperature structure in 
the intervening decade. The disk-averaged nature of the 
AKARI spectra prevent detection of spatial variations in 
the temperature field, though we expect these conclusions 
to be weighted towards the warm south polar region of 
Neptune where th e hydrocarbon emissions are strongest 
( Qrton et al. 20071 ) . The CH4 abundance is larger than 
that expected if the tropopause at 56 K is efficiently cold- 
trapping methane in the troposphere, but consistent with 
the hypothesis that CH4 lea king through the w arm south 
polar tropopause (62-66 K, lOrton et al. 20071 ) is globally 
redistributed by stratospheric motion. 

The uncertainties on the retrieved T(p) profile are large 
because of the strong anticorrelation of temperature and 
CH4 abundance. However, this was subsequently used to 
derive (a) the ratio of D/H in CH4 of 3.0 ± 1.0 x 10"'' from 
the CH3D emission at 8.6 fim; (b) the stratospheric mole 
fraction of ethane of (8.5±2.1) x 10"'' at 0.3 mbar from 12.2- 
/im; and (c) the mole fraction of ethylene of 5.0^2^ x 10~^ 
at 2.8 /ibar from the 10.5-/im emission band. The D/H ra- 
tio supports the conclusion that Neptune was enriched in 
deuterium relative to the other gas giants at th e time of 
the planets formation (|Feuchtgruber et al. 1999f ). The ob- 
served ratios of the hydrocarbon mole fractions will pro- 
vide important inputs to constrain photochemical model- 
ing for Neptune's stratosphere, and explain the observed 
differences between the giant planets. 

At shorter wavelengths, Neptune flux densi- 
ties in the 3.5-4.1 /zm range (indicative of the 
production of stratospher ic hydrocarbon hazes, 
iFeuchtgruber fc Encrenaz 2003h show considerable varia- 
tion with time. Comparison with previous studies confirms 
that the flux density is highly correlated with solar 
variability (determined from the solar sunspot cycle), with 
the low flux measured by AKARI the direct result of the 
present solar minimum. AKARI spectra demonstrate the 
first evidence for Neptune's fiuorescent emission of CO 
at 4.7 /im, comparable with similar emission observed 
recently on Uranus (Encre naz et al. 2004). The fluorescent 
emission can only be reproduced using a vertical CO 
distribution simil ar to that derived fr om sub-millimeter 
wavelengths by iHesman et al. f2007T ). confirming the 
necessity for two sources of CO in Neptune's atmosphere, 
one from convective mixing with the deep troposphere, 
and another exogenic source for the stratosphere. 
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Table 1. Akari/IRC Neptune Data 



Module 


UT Time 


Subsolar Longitude 


Wavelength 


Resolution 


NG 


1756-1826 


165 


2.5-5.0 fim 


135 at 3.6 /xm 


NP 


0302-0332 


189 


1.7-5.5 


22 at 3.5 fim 


SGI 


0302-0332 


189 


5.5-8.3 pim 


47 at 6.6 ^im 


SG2 


0302-0332 


189 


7.4-13 ^lm 


34 at 10.6 



Table 2. Correlation between L-band flux density and sunspot number (SSN) . L-b and flux densities wer e es- 
timated from the spectra sho wn in the relevant publications [1] iPilcher (19771 ). [2] lEncrenaz^ et al. ('2nnnal) . [3] 
iFeuchtgruber fc Encrenaz (^2003f ) and [4] This work. Monthly averages and standard deviation for SSN were taken 
from th e International Sunspot Numbers, compiled by the Solar Influences Data Analysis Center in Belgium 
( http: / / solarscience.msfc.nasa.gov/greenwch/spot-num.txt, ). 



Date 


L-Band Flux Density 
mjy 


SSN 


Reference 


May 1975 (Kitt Peak) 


< 10 


9 ± 11.1 


[IJ 


May 1997 (ISOPHOT) 


7 


18.5 ± 13.4 


[2] 


Aug 2002 (VLT) 


20 


116.4 ± 36.3 


[3] 


May 2007 (AKARI) 


5.5 


11.7 ± 6.8 


[4] 
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Fig. 1. AKARI/IRC Spectrum of Neptune from 13 May 2007. The upper panel plots the radiance per wavenumber, 
which emphasizes the near-IR portion of the spectrum from the NP (circles) and NG (triangles) modes. The lower panel 
plots the radiance per wavelength, which emphasizes the mid-IR spectrum from the SGI (crosses) and SG2 (diamonds) 
modules. SGI data between 5.5-5.8 /im and NP data shortward of 5.0 /im are omitted for clarity. 
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Fig. 2. Contributi on function for the SGI and SG2 spectral range, calculated using the reference atmosphere of 
iMoses et al. (20051 ) . The darkest colours indicate the peaks of the contribution function. Note in particular the multi-lobed 
nature of the CH4 contribution function at 7.7 /im, used to derive the stratospheric temperature structure. 




Temperature (K) Temperature (K) 

Fig. 3. Family of models (a-d described in the main text) for the derivation of Neptune's quasi -isothermal stratosp heric 
temperatures. The T{p) profile for pressures greater than 1 mbar was set to the profile used bv lMoses et al. (2005[) . and 
the T{p) in the 1-1000 /ibar region was varied simultaneously with the stratospheric CH4 abundance to fit the 7.7 /im 
emission band. 
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(a) Variation of Methane Abundance 
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(c) Variation of Temperature 



(b) Trade-Off Analysis 
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Fig. 4. Demonstrating the trade off analysis for fitting temperatures and methane with the family of models in Fig. [3l 
Panels (a) and (c) show the competing effects of varying CH4 and T{p), and panel (b) shows how the best fitting model 
and error range was estimated from a 2-dimensional analysis. 




Fig. 5. Best-fitting temperature profile (solid line) compared to the results of a selection of Neptune analyses. The 
grey shading indicated the range of possible solutions in the stratosphere given that the CH4 abundance is highly 
anticorrelated with temperature (see main text). Nevertheless, the tradeoff analysis in Fig. 2] provides stratospheric 
temperatures consistent with previous analyses from ISO and Voyager. 
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Fig. 6. Varying the abundance of CH3D has a small effect on the spectrum at 1160 cm~^. Although the spectral lines of 
CH3D are not observed, we can broadly fit the shape of the smooth spectrum to estimate the D/H ratio in stratospheric 
CH4. The steps are in increments of 2 x 10^^ on the D/H ratio. 




780 800 820 840 860 880 900 
Wavenumbers (cm ') 

Fig. 7. Plotting the variation of the modelled radiance compared to the AKARI data for the C2H6 e mission feature. Th e 
thick solid line is the model with the best-fit T{p) from Fig.[5]and 0.8 times the modelled abundance of (jMoses et al. 2005D . 
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Fig. 8. Plotting the variation of the modehed radiance compared to the AKARI data for the C2H4 emission feature at 
10.5 fim. The signal was extremely small between the much larger ethane and methane emissi on features, but the thick 
solid line is the model which best-fits the morphology of the spectrum, scaling the model of iMoses et al. (2005[ ) by a 
factor of 1.4. 
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Fig. 9. Fluorescence of CO compared to the near-IR prism spectrum from AKARI. The thick solid line is the best fitting 
profile, with CO present in both the troposphere and stratosphere (profile 2). The dashed line is the modeled emission 
with stratospheric CO only, indicating the requirement for tropospheric CO. 
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Fig. 10. Comparison of CH4 (circle, compared to the solid line), C2H6 (diamond, compa red to the dotted lin e) and C2H4 
(triangle, compared to the dashed line) results with the predictions of Model A from iMoses et al. f2005f ). The model 
predictions for the vertical distributions, which were shown to be consistent with Voyager/UVS results in Fig. 32 of 
Moses et al. f200"5l) . were scaled to fit the AKARI spectra. The do tted box shows the range of CH4 mole fractions from 
Yelle et al. the dashed box shows the range of values from lBishop et al. (1992t) . 



